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Abstract 

We consider the homogeneous cosmological radion, which we define 
as the interbrane distance in a two-brane and Z2 symmetrical config- 
uration. In a coordinate system where one of the brane is at rest, 
the junction conditions for the second (moving) brane give directly 
the (non-linear) equations of motion for the radion. We analyse the 
radion fluctuations and solve the non-linear dynamics in some simple 
cases of interest. 
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1 Introduction 



Motivated by recent developments in string theory, the idea of confin- 
ing matter in a three-brane embedded in a higher dimensional space- 
time has stimulated lately a lot of activity in several fields of research. 
In cosmology, it has been shown that, quite generically, the matter 
confinement in a brane entails a non-standard evolution of the brane 
geometry |jl|. In other words, the Friedmann equations governing the 
expansion of our 'apparent' universe are modified. This could lead to 
trouble because nucleosynthesis would give very different light element 
abundances with a modified scale factor evolution. Therefore, a strong 
constraint on brane models is the recovery of standard cosmological 
evolution at least since nucleosynthesis. A nice solution to this prob- 
lem has been provided by applying the idea of the Randall- Sundrum 
model m that a tension in the brane can be compensated by a neg- 
ative cosmological constant in the bulk, as far as the internal brane 
expansion is concerned. Applied to cosmology, this idea leads to the 
recovery of the usual Friedmann law when the ordinary cosmological 
energy density is much less than the brane tension Q. By contrast, 
at earlier times, when the energy density is higher than the tension 
one finds a non-standard evolution in the brane [^. 

Another solution that has been suggested in the case of finite size 
extra-dimension (see e.g. Q for the case of infinite extra-dimension) 
is to consider a stabilizing potential for the so-called radion ||^ ^. The 
precise definition of the radion is the source of some confusion in the 
literature. It is usually described as the metric component along the 
extra-dimension(s), but as emphasized in it can be reinterpreted 
as the relative distance between the two branes. Stabilization of the 
radion can be explicitly realized [10| by adding some matter in the 
bulk, a scalar field for simplicity, so that the relative distance between 
the two branes will be stabilized. This mechanism has been applied 
[ p^ , PH , p^ , p^ ] to the two-brane model of Randall and Sundrum 0, 
in order to recover gravity consistent with the observations on the 
negative tension brane (without stabilization, one gets a Brans-Dicke 
gravity with a negative parameter ||l^ ) . In cosmology, this mechanism 
has been discussed in [15] and |16| but the claim that this would lead 
to the usual Friedmann equations has been challenged in |17]. Let 
us mention that the radion has also been discussed in the context of 
holography |jTM|. 
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The purpose of the present work will not be to study the cosmo- 
logical consequences of a stabilizing mechanism for the radion, task 
which will be left for future work, but more modestly to characterize 
the radion in cosmology (see also a brief discussion in to write 

down the equations of motion that govern its dynamics, which will be 
analysed in some simple examples. 



Previous treatments of the radion in brane systems 12, 20, 21 
start from the beginning with a perturbative approach: the radion 
is thus hidden among the metric fluctuations, which have been the 
subject of active research lately in cosmology (see |2^ and references 
therein). Our approach is quite different in the sense that the radion 
dynamics appears directly from the junction conditions. It also gives 
us the full non-linear dynamics of the homogeneous radion (i.e. ig- 
noring the ordinary spatial dependence of the radion). We are thus 
able to show that the dynamics of the radion is more subtle than that 
of a scalar field coupled to the trace of the matter energy-momentum 
tensor. 

Our plan is the following. We begin in section 2 by defining the 
coordinate system in which we describe the two-brane configuration. 
We then focus on the internal geometries inside the two branes in sec- 
tion 3. Section 4 is the central part of this paper, where the equations 
of motion of the radion are derived. Starting in section 5, we apply 
our general results to the case of an empty anti-de Sitter (AdS^) bulk 
spacetime. The last two sections are devoted to the analysis of the 
radion evolution in this geometry, first using a perturbative approach 
(section 6), then in the full non-linear case (section 7). 



2 The two-brane cosmological system 

The purpose of this section will be to present in a geometrical way 
the space-time configuration we wish to study. The spacetime will 
be supposed to be five-dimensional and to satisfy three-dimensional 
spatial homogeneity and isotropy in order to obtain (unperturbed) 
cosmological configurations. Two spatially homogeneous and isotropic 
three-branes live in this spacetime. The fifth dimension is periodic 
and we assume a mirror (orbifold) symmetry across each of the branes. 
The space between the two branes represents half of the periodic space 
along the fifth dimension, the other half being the mirror image with 
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respect to any of the two branes. 

In any coordinate system, even adapted to the three-dimensional 
homogeneity and isotropy, the two branes would be in general in mo- 
tion along the fifth dimension. It turns out that it is possible, and 
convenient, to define the fifth dimensional coordinate y as part of a 
Gaussian normal coordinate system with respect to one of the branes, 
chosen arbitrarily and called the reference brane. The reference brane 
is at rest in the fifth dimension, in y = say, and the spacetime metric 
in this coordinate system (which is not unique) is of the form 

ds"^ = gAsdx^dx^ = —n{t, y)'^dt^ + a{t, y)'^^ijdx'^dx^ + dy^, (1) 

where the are coordinates for the maximally symmetric three- 
dimensional submanifolds, t is a time coordinate. The metric com- 
ponents are necessarily of this form because of the symmetries, but 
their explicit dependence of t and y can be obtained only by solving the 
five-dimensional Einstein's equations (with the appropriate boundary 
conditions). 

Gab = Rab - ^Rqab = k^Tab, (2) 

where Rab is the Ricci tensor, R its trace, and Tab the energy- 
momentum tensor. The coupling of the matter to the geometry, k^, 
can be written as the inverse of the third power of the five dimensional 
reduced Planck mass. 

If the reference brane is now at rest, the second brane will be in 
general moving in the same coordinate system. Its position at any 
time t will be denoted 

y = n{t), (3) 

and we call this function TZ{t), which represents the relative distance 
between the two branes, the cosmological radion. Note that the role 
of the reference brane and of the second brane can of course be ex- 
changed. Finally we wish to stress that here, unlike most of the liter- 
ature, our radion is defined in a non-perturbative way. We will show 
later that our definition coincides, in the linearized approach, with the 
traditional definition of the radion. 

3 Induced geometries in the two branes 

The induced metric in the reference brane immediately follows from 
the spacetime metric (^ and can be written in the usual Robertson- 
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Walker form, 

dsfnduced = -dt^ + aldx^, (4) 

where the time coordinate has been redefined such that n(t, 0) = 1, 
and ao{t) = a(t, y = 0). t is thus the proper time associated with the 
reference brane. 

To get the induced geometry for the second brane, one must take 
into account its displacement in the fifth dimension, and one finds 
from (0) the induced line element 



"''^induced 



n2(^,7^(^)) -Tt^ dt'^ + a'^{t,n{t))dx^, (5) 



where a dot will always stand for a (possibly partial) derivative with 
respect to t. This can be rewritten in the form 



dsinduced = -dr + a2{T)dx , (6) 

where r is the second brane cosmic time (i.e. the proper time for 
a comoving observer) which is related to the reference brane cosmic 
time t by 



dT = n{t,nit))Jl--j^pj^ dt = n2^'Ut. (7) 

And the scale factor in the second brane is simply 

a2{t)=a{t,n{t)), (8) 

or 0^2 ("T") if one expresses t in terms of r. 

It is also useful to introduce the expansion rate for the second 
brane, which can be defined as 

H2it)^^ = (^ + ^n) . (9) 

However, H2 does not coincide with the standard definition of the 
Hubble parameter for a cosmological observer in the second brane be- 
cause it is defined with respect to the time t, which is the cosmic time 
in the reference brane but not in the second brane. The traditional 
Hubble parameter for the second brane corresponds to the definition 

^^^^^ ^ = :lh2 (10) 

a2 n2 
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Finally, one can introduce the unit velocity vector corresponding 
to the second brane, whose components read 



A r dt dy . 1 , dTZ , 

= 0, 0, 0, ^} = n2-S{l, 0, 0, 0, — }. (11) 

At this point, we can easily infer the cosmological evolution in the 
second brane from the knowledge of the cosmological evolution in the 
first brane and the position of the second brane with respect to the 
first. What will be studied next is the connection with the matter 
content of the second brane, which has been ignored up to this point. 



4 Junction conditions 

The matter content of the branes is directly related to the jump of the 
extrinsic curvature tensor across the brane. This relation has been 
established explicitly in the case of a brane at rest with respect to 
the coordinate system (||) in our previous works [||, ^. We generalize 
our result to include the case of a brane moving with respect to the 
coordinate system. 

We will now consider explicitly the second brane, but all the follow- 
ing expressions apply as well to the reference brane by taking TZ = 0. 
The extrinsic curvature tensor on the brane is defined by the expres- 
sion 

Kab = h^VcriB, (12) 
where is a unit vector field normal to the brane worlsheet and 

hAB = QAB - riATlB (13) 

is the induced metric on the brane. 

It is easy to compute the components of the unit normal vector 



from the components of the brane velocity (11), by using the fact that 
n"^ must satisfy the two following conditions: 

QABn'^n^ = 1, gABTi^u^ = 0, (14) 

and we get 

= { ^=^'0'0'0' ; V J . (15) 
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Then, a straightforward calculation consisting in the substitution of 



the components (|l5| ) into the definition (12) yields the following (non 
zero) components for the extrinsic curvature tensor: 

< = : ^ (16) 



Kl = TZKl Kl = -Kl/n\ (17) 
K] - ^f^ + Sk, (18) 




{1Z + nn' - 2^7t2 - ^7^ 
"572 



= ^ . ^..n " ^ (19) 



where all coefficients take their value on the brane, i.e. at t and 

y = n{t). 

Let us now introduce the energy-momentum tensor on the second 
brane, which, because of the symmetries of the setup, is necessarily of 
the perfect fiuid form 

SaI = (P2 + P2)UAUB + P2hAB. (20) 

We also define Sab as 

Sab = Sab - -^ShAB- (21) 



The junction conditions for a hypersurface in a five-dimensional 
world read [||] 

[K^.] = -z^'^^., (22) 



where the brackets stand for the jump across the brane, and 
e^e^KAB (and a similar expression for Sf^^) where the four vectors 
{fj, = 0, 1, 2, 3) form a basis of the vector space tangent to the brane 
world volume. For the second brane, this gives explicitly 

[K^,,]^ = K^,{t,TZ{t)+) - K^,{t,n{t)-) = -2K^,{t,n{t)), (23) 

where the last equality follows from the identification of the point 
y = TZ{t) with the point y = —TZ^t) and from the Z2 symmetry across 
the second brane. Note that, whereas all expressions (15-22|) apply also 
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to the reference brane (by taking 7^ = 0), the last expression differs 
by a sign for the two branes: one would find [-f^/xi/Jo = '^Kf^y{t,y = 0) 
for the reference brane. 

It is easy to compute the components of S^, using the velocity 
components (|Tl|) and one finds 



Si = J_P1±^^^ (24) 



St = ^Pl±^^2^ Si = -''-P^±^^'n (25) 
1 

Si ^ ^^,4. (27) 
3 n"^ 

Equating this with the extrinsic curvature tensor, according to (22), 
all equations reduce to only two equations, which read 

2 / -772X^/2 



S] = ^925] (26) 




nil k\ 7^■ 

!— 2 = "IT 2P2 + 3P2 1-- 




a' alZ ■ - I ' - 1 / N 

- + -— = -^2 1-— , 29 
a an'^ " 

where, once more, all the metric coefficients take their value on the 
brane. The above system describes the full non-linear dynamics of 
the homogeneous cosmological radion. The non-linearity is manifest 
in the terms quadratic in TZ but is also hidden in the fact that all the 
metric coefficients depend on IZ. 

It is interesting to note that the left hand side of equation (|28| ) 
is simply the (covariant) acceleration in the fifth direction of a par- 
ticle with the velocity (|lT|). This can be understood by noticing 
that u'^Kabu^ = —riBU^Vcu^, as a consequence of (12) and (p^). 
u^Kabu^ : which is also —Kq, therefore vanishes for a geodesic mo- 
tion. The right hand side of equation (^) will thus induce a deviation 
from geodesic motion for a "comoving" observer on the second brane 
(at fixed coordinate x*). 

The second order equation (|28| ) can in fact be derived from the first 
order equation (p9|). Indeed, if one takes the time derivative of ( [29| ) 



and adds it to H2 times (29), one gets, after an overal multiplication 
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by 7 



2, the expression 




(30) 



where we have used 



P2 = -3i?2 (P2 + P2) 



(31) 



which expresses the energy conservation on the second brane (assum- 
ing there is no inflow/outflow of energy from/to the bulk). A bulk 
energy-momentum tensor, satisfying the three-dimensional symme- 
tries of homogeneity and isotropy, can be written in the form 



where v is a unit time-like vector (whose ordinary spatial components 
are zero) describing the fluid velocity, Pt represents the transverse 

pressure (along the fifth dimension) and Pb represents the pressure in 

(2) 

the three ordinary dimensions; is the projection tensor on the two- 

dimensional (t, y) subspace and is the projection tensor on the 
three-dimensional ordinary space. Because of the Z2 symmetry, the 
bulk fluid cannot go "across" the branes, and we must thus impose the 
further condition on the bulk fluid that the vector field coincides, 
on the branes, with the brane velocity li"^ Q . It is then not difficult to 
check explicitly that the second and third terms on the left hand side 
of (|30|), expressed in terms of Tab by use of Einstein's equations (Q), 
just cancel. Equation ([2^) is thus derived, provided H2 is not zero. 

Equation ( p9|) can be solved algebraically to express TZ in terms of 
the other quantities. One finds 




(32) 




Simultaneously + vavb — nAns- 
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If a! I a and p2 are of the opposite sign, then (^) has one solution if 
|d/(an)| > |a'/a|, no solution otherwise. If a! ja and are of the same 
sign, then p9| ) has two solutions when \aj{an)\ < \a' /a\, and only one 
solution otherwise. In both cases, when there is a unique solution, it 
corresponds to the root in ( |3^ ) with the sign being the same as that 
of d/a. 

To conclude this section, let us mention that the equations of mo- 
tion for the radion can be rewritten in terms of the cosmic time on the 
second brane, r, instead of the cosmic time on the reference brane t : 

-^{2p2 + W2)^l + n% (34) 

tp^(i + niy"\ (35) 

The second equation has been obtained by |2^] in the particular case 
where the reference brane is Minkowski (hence d/a = and p2 is a 
constant. 



n^, + - f 1 + n\ 

n V 



',rr 



+ 



a 



an 



5 Radion dynamics in an AdS^ bulk 
spacetime 

In order to explore the dynamical evolution of the radion, we now 
need to be more specific and to resort to explicit expressions for the 
spacetime metric (|l]). We will thus choose to work with an AdS^ 
bulk spacetime, for which it is possible to recover (approximately) 
at late times, standard cosmology, and which corresponds to the 
cosmological, i.e. non static, generalization of the Randall- Sundrum 
model 1^]. 

Keeping aside the two branes, we assume that the bulk spacetime 
is empty, except for a negative cosmological constant A. It is useful to 
define a mass scale p associated with A, according to the expression 

p = y^A/e, (36) 

{p~^ is the AdS^ curvature radius) and an energy density scale a, 
defined by 
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Note that, in the Randah-Sundrum model Q, the energy density in 
the brane is supposed to be exactly this a. In the present work, we 
consider any value for the energy densities in the two branes. 

We will also assume that the five-dimensional Weyl tensor is zero. 
Our global spacetime is therefore made of two pieces of AdS^ separated 
by the two branes. Assuming the usual Z2 symmetry with respect to 
each brane, the two AdS^ pieces are then simply mirror symmetric. 

In Q, we have solved explicitly the five-dimensional Einstein's 
equations 

Gab - J^QAB = i^^Tab, (38) 

in the case where the only matter in spacetime is the one confined in 
a brane. We obtained the following form for the component of @): 

a{t,y) = ao(t) (cosh/iy - 7yosinh/x|y|) (39) 
n{t,y) = cosh/xy — r^o sinh/x|i/|, (40) 

where rjo and fjo are two functions of time, related to the energy density 
in the reference brane by 

Vo{t) = , Vo{i) = % + TT' v^^) 

a Ho 

where Hq = do/ao is the Hubble parameter. For future use we also 
define r/2 as 

V2{t) = (42) 
a 

where P2{t) is the energy density of the second brane. It must be 
stressed that, whereas the functional form of the metric components 
( |39| , pC| ) is obtained by solving the Einstein equations in the bulk, the 
relationship between the coefficients and the reference brane matter 



given in (41) comes from the junction conditions |l[ across the refer- 
ence brane. 

Finally, the evolution of the scale factor ao(^) and of the energy 
density po{t) is determined by solving the coupled system consisting 
of the generalized Friedmann equation 

where k = —1,0,1 according to the three-dimensional spatial cur- 
vature (the term C/a^ does not appear because we have assumed a 
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vanishing Weyl tensor) and the usual energy conservation law 



po = -3Ho {po + Po) . (44) 
A useful identity is an integral of Einstein's equations 



77^^^^+ - (45) 



na) 6 V fl / « 

which was obtained in |p. The generalized Friedmann equation ( |4^ ) 
on the reference brane is simply the application of this identity (with 
the junction conditions) in y = 0. 

The fully non-linear evolution of the radion, coupled to the evolu- 
tion of the matter on the two branes, is thus given by the first-order 



system of differential equations consisting of (|43|), (|44D , (|3lD and (|33|), 
which can be rewritten in the simpler form 

where 7^2 is the Hubble parameter of the second brane defined in ( p!o| ) . 

One can also check, starting from the expression (P) in terms of 
the metric component ( p9| ) and using ( [29| ) and (|^), all terms in 7t 
cancel and one recovers the generalized Friedmann equation on the 
second brane, that is 



ni = -^^^-p.'-K = P?{r,l-\)-K. (47) 
6 36 a2 ^ '0-2 



It has to be so because the generalized Friedmann equations comes 
from a purely local analysis, and the role of the two branes can be 
exchanged. 



6 Radion fluctuations 

In this section, we study perturbatively the dynamics of the radion 
fluctuations. We first define the background solution as the case where 
the radion is frozen, i.e. does not depend on time. This means that 
the second brane is at rest with respect to the first brane, i.e. 

y2(t) = % (48) 
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where is a constant. One can then find exphcit solutions of the 
Einstein's equations for the whole two-brane system. However, to 
ensure the continuity of the metric along the compact fifth dimension, 
one finds that the matter content in the second brane is necessarily 
related to the matter content of the first brane. This result was already 
noticed for a vanishing bulk cosmological constant Q. We generalize 
it here to the case of a non-zero bulk cosmological constant. 

If TZ is constant, then the system ( p8[]2^ ) simplifies considerably to 
yield 



. — {2p2 + 3P2) (49) 
n J o 6 



and 

/a'\ , , 

Substituting the explicit expressions (^9|) and (|^) for a and n, which 
depend on the energy content on the first brane, one can express the 
energy density and pressure of the second brane as functions of the 
energy density and pressure on the first brane and of the position TZ. 
One finds 

2-(.\ r sinh /i7^ - r/o cosh /i7^ 

K p2{t) = 6p — ^ ■ . ^ , (51) 

cosh fj,l<. — riQ smh ///<. 

and 

2- sinh/x7^-f?ocosh^7^ 2 2- 

K P2{t) = -2/i — ^ _ - -K p2. (52) 

cosh fil-c — ijQ smh pl-c 6 

Note that if one considers an equation of state of the form Pq = 
wqPo with wq constant for the matter on the reference brane (which 
includes the standard cases of radiation, non-relativistic matter and 
a cosmological constant), the equation of state of the matter on the 
second brane, which is completely determined here as a function of 
the relative distance between the two branes, will correspond to an 
equation of state P2 = W2P2 where, in general, W2 is time dependent. 
Indeed r]o and 

% = -(2 + 3wo)% (53) 

are time dependent in general. 

One exception is when the reference brane contains only a cosmo- 
logical constant, i.e. wq = —1, then one also gets 102 = —1 on the 
second brane. There is in this case an interesting correspondence be- 
tween the cosmological constants on the two branes. Introducing the 
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effective cosmological constant from the brane point of view as 



A 

^=6+36^ 



1), 



(54) 



we have 



An 



v2 • 



> 1) 



(55) 
M4 



(cosh ^TZ — r/o sinh ^iR) 

Thus, the two branes are simultaneously dS^ (A > 0, \rj 
(A = 0, |?7| = 1) or AdSi (A < 0, |r/| < 1) ||. 

Another interesting case is when the reference brane undergoes 
ordinary cosmology, which occurs when t/q is very small. One can 
then neglect 770 as well as 7?o in the expressions (^52) and the matter 
on the second brane reduces effectively to a cosmological constant. 

In order to obtain the equation of motion for the radion fluctu- 
ations, we linearize the exact equations of motion ( p8[]29| ) about the 
equilibrium configuration we have just defined. The linearized system 
reads 



5n 



— + 5 ( - I ^ 

h5TZ 
+ — 



6\- I + J 



-— (25p2+35P2) (56) 
o 

y5p2, (57) 



where we recall that all metric coefficients take their value on the 



second brane. Furthermore, one can use the explicit expressions (39 
1^) to find 



5 - 



m^5TZ 



ml5TZ 



36 

4 



(2/32 + 3P2) 



2 ^ - 2 

^ "36^^ 



5n, 



(58) 
(59) 



where we have introduced two effective squared masses and m^. 
Note that both m„ and rua vanish in the particular case where the 
Randall- Sundrum condition between the tension and the bulk cosmo- 
logical constant is satisfied on the second brane. 

If one adds (^) with three times (|57|), one obtains the equation of 
motion 



^ + (;3«-^)^ + (m^ + 3m5^7^ = 4^T2, 
\ a n J n-^ \ ' 6 



(60) 
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where 5T2 = —Sp2 + 3(^P2 is the hnear perturbation of the trace of the 
matter energy- momentum tensor. One thus recognizes the equation of 
motion for a scalar field coupled to the (perturbed) trace of the matter 
energy-momentum tensor. The friction coefficient is the usual one 
(rewriting in terms of r would give the natural 37^2 friction coefficient, 
as can be seen when extracting this equation directly from (^), (pSD), 
but the mass term is an effective mass that depends both on the bulk 
cosmological constant and on the background matter content of the 
second brane. 

This effective mass is time-dependent except in the case where the 
second brane matter has for equation of state P2 = —p2- The two 
squared masses are then identical and equal to 



m„ 



-A2 



(61) 



where A2 is the effective cosmological constant from the brane point 
of view defined in (H). One thus recovers exactly previous results 



2C, 21]. An immediate consequence of this is that when the first 



brane is dS^^ the small perturbations around the equilibrium position 
of the second brane (then also dS^) are unstable, whereas when the 
first brane is AdS^ (the second brane is then also AdS^i) they are 
stable. We will see in section |^ that this perturbative analysis matches 
with the exact non linear evolution. 

The evolution of the radion fluctuations is entangled with the mat- 
ter fluctuations. It is however possible when the matter fluctuations 
are of the adiabatic type, i.e. 



= Cs2Sp2, 



(62) 



to derive an equation for only the radion fluctuations by considering 
the appropriate linear combination of ( |5^ ) and (^), 



sn 



+ 



(2 + 34) 



n 



5n + 



ml +{2 + 30^2) "la 



5TZ = 0. 



(63) 

The evolution of the matter fluctuations is then related to the radion 
fluctuation by (57), i.e. 



—6p2 = 2" + ^aOT^- 

6 a 



(64) 
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Assuming an equation of state P2 = W2P2 (and adiabatic fluctua- 
tions, i.e. c^2 = '^^2), the equation of motion (|6^ ) takes the form, after 
substitution of and by their explicit expressions, 

Sn^rr + (2 + 3W2)n2dnr + 3(1 + W2) " |^/'"2^(2 + 3'u;2) ^ STZ = 0, 

(65) 

where we have used here the second brane cosmic time r. One can 
solve this equation for de Sitter {w2 = —1) on the two branes and one 
gets 

6n = 57^,e^^(---») - ^^6p2^ (l - e^^(---')) , (66) 

where 6TZi is the initial radion fluctuation and 5p2i is the initial density 
fluctuation (at time Tj). For de Sitter brane cosmology, we recover the 
instability mentioned above. 

Note that in the case of the radion fluctuations in the strict Randall- 
Sundrum model, studied in Q, the friction coefficient in the above 
equation vanishes because the background is static, and both the ef- 
fective squared masses and also vanish because of the RS con- 
dition relating the bulk cosmological constant and the brane tensions. 
One thus ends up with the equation of motion for a free scalar field. 



7 Non linear evolution of the radion 

In this section, we explore the non-perturbative aspects of the radion 
evolution in some simple cases, leaving a more systematic analysis for 
future work. We are going to consider the cases where the matter 
on the reference brane behaves like a cosmological constant. In these 
simple cases, the reference brane matter does not evolve, i.e. rjo = 0, 
and the metric components are separable and can be written in the 
form 

n{t,y) = n{y), a{t,y) = ao{t)n{y). (67) 

The explicit form of n{y) depends on the sign of t]q — 1. If rjQ > 1, 
the brane is d5'4; if ?7o < 1, the brane is AdS^; and the case r]o = 1 
corresponds to the Randall-Sundrum model [^. We will then fur- 
ther specialize to the cases where the brane matter of the second also 
behaves as a cosmological constant, implying 7)2 = 0. 
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7.1 dS^ brane 

For a (i5'4 reference brane, the dependence of the metric components 
on the fifth coordinate is simply (see e.g. p^]) 



'^(y) = V^o - lsinh/i(2/h - |y|), (68) 

where ^ 

yh = - arg coth r]o (69) 

Here we suppose that the second brane is always located at 7^ < yh- 
Otherwise, there would be a horizon at 7?. = yh- 



Equation (33) can then be rewritten in terms of dimensionless 



quantities in the form 

2 ho cosh x±ri2 sinhx-y/ r^^/Cg sinh^ x — /Cq cosh^ x + 

= -/Cosinhx ————-^^ 

(70) 

where we have introduced the dimensionless distance of the second 
brane to the horizon x and the dimensionless time s 

X = ii{yh - 7^), s = /it. (71) 

The dimensionless Hubble parameter /iq is defined by 

ho^-^ = ^, (72) 

and /Co is defined by 

/Co = ^\vi-l\ (73) 

In general equation ( |70| ) is non separable since both 7/2 and ho are 
function of time. It gets however simplified when one assumes that the 
energy momentum tensor of the second brane is that of a cosmological 
constant (implying 7)2 = 0), and if one assumes moreover that /c = (in 
order to account for a slicing of dS^ by flat 3-spaces) so that ho = /Cq 
(as can be seen from (^3|) ). In this case equation ([70|) reads 



cosh X ± r/2 a/t/I ~ 1 sinh^ x 

x., = -\/r]^-lsmhx — — . ^ ' (^4) 

1 + smh X 
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it becomes separable and can be solved analytically. Since the matter 
energy densities are frozen on both branes, the evolution of the radion 
is fully embodied in this unique first order equation. The evolution 
of the radion is thus completely determined by the knowledge of the 
initial position of the second brane, which we call TZi. If the initial 
position TZi is located at the equilibrium position 



then the brane does not move. However, if the brane is initially away 
from the equilibrium position, it moves. According to the perturbative 
analysis, a brane slightly off the equilibrium position moves further 
away since we found that the radion is unstable. 

Here, we can determine the non-linear evolution. We have solved 
numerically the evolution of the radion in the cases where the brane 
is initially on the left {TZi < TZ) and on the right (TZi > TZ) of the 
equilibrium position. We observe the expected unstable behaviour. 
In the first case, illustrated in Fig. |^, the brane moves towards the 
reference brane until it collides with it (numerically it goes across it). 
In the second case, illustrated in Fig. |2|, the brane moves further away 
but its motion freezes asymptotically at the horizon, as seen in the 
reference brane proper time. However, as can be seen from (|^) and 
([74[), this takes only a finite amount of the second brane proper time, 
as is typical for the crossing of a (Rindler or black hole) horizon. 

One may obtain a simple analytical expression which illustrates 
this behaviour in the case of a empty second brane : ?72 = 0. The 
solution is 



where Xi is the initial position at time Sj. This is also the solution of 
( [74[ ) for non zero r/2 when the second brane approaches the horizon 
(i.e. for X small). 



For an AdS4 reference brane one finds for the metric component 



X = fi{yii - TZ) 



arg coth r]2 



(75) 




(76) 



7.2 AdS^ brane 




(77) 
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where ^ 

yc = - argtanh?7o. (78) 

In this case the warp factor n has a minimum in y = yc enabling to 
have a configuration with two static branes of positive tensions [25| 



(see also |21, 27|). In any case, equation (p3) reads 



— /lo sinhx ± r/2 cosh x^t/I/Cq cosh x — /Cg sinh x + 
Kn cosh X 



/iq + ??2^o cosh^ X 
where x is defined by 

x = ^{TZ-y,), (80) 
and s and /Co are defined respectively as in (0) and dfl). The 



(79) 



scale factor is given straighforwardly by solving (|4^); one finds ag = 
sin(/Cos)//x/Co, so that in this FLRW parametrization of AdS^ there 
are coordinate singularities in ICqs = and ICqs = vr. The dimen- 
sionless Hubble parameter is then /iq = /Cocotan(/Cos). The equation 



of motion (7£) remains non separable even in the simple case where 
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= 0, with r/2 7^ 0. We have solved numerically for the evolution 
of the radion in the cases where the brane (with r/2 being constant 
and positive) is initially (at time such that ICoSi < vr/2) on the 
left {xi < x) and on the right {xi > x) of the equilibrium position x 
defined by 

X = ^{TZ — Uc) = argtanhr/2. (81) 

For small enough Si, the only acceptable solution is given by taking 
the + sign in front of the square root in equation (|7^) , as can be seen 
using the arguments given at the end of section ^ 

In both cases {xi < x and Xi > x), illustrated in figures ^ and ^ 
the brane first moves back toward its equilibrium position, and then 
reaches a point (at time Sc close to 7r/2/Co) where the square rootQ in 



equation (^9|) vanishes. The solution obtained from ( |79| ) with the + 
sign in front of the square root must then be matched (after time Sc) 
to a solution of (|79|) with a — sign is taken in front of the square root 
(the two solutions have the same tangent in Sc). Figures § and § are 
obtained with this matching. 

Equation (|7^ becomes separable, and can be easily integrated, in 
the limiting case where the second brane is empty i.e. r]2 = 0. In this 
case, one finds: 

cos (JCqs) 



X = argtanh 



tanh X- 



cos (/CoSj) 



(82) 



7.3 M4 brane 

When the reference brane is Minkowskian (t/q being 1, and a = n = 
exp(— /i|y|)) , equation ( |33| ) is given by 

x,s = ±^^V2-h (83) 

where x is given here by 

x = fin (84) 

and s is defined as above. The second brane equilibrium condition 
is simply given by the Randall- Sundrum [^, |^ matching condition 
expressed here as r/2 = —1- For this particular r]2, one sees easily 

^which is also the determinant of the quadratic equation obtained from equation 
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Figure 3: Evolution of the distance x — x between the brane position and 
the equihbrium position as a function of /Cqs/tt, in the case of AdS^ for an 
initial position Xi < x. 
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Figure 4: Evolution of the distance x — x between the brane position and 
the equihbrium position as a function of /Cqs/tt, in the case of AdS^ for an 
initial position Xi > x. 
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from equation ( |83| ) that the second brane must be at rest with respect 
to the first one, but that the equihbrium position x can be chosen 
arbitrarily. However a small departure from the equilibrium condition 
r/2 = — 1 results in either a colliding of the two branes or a runaway 
behaviour. 
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